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ABSTRACT _ 

This  paper  describes  the  design  and  the  simulation  of  a  non-linear  controller 
for  two-mass  system  using  induction  motor  basing  on  the  backstepping 
method.  The  aim  is  to  control  the  speed  actual  value  of  load  motor  matching 
with  the  speed  reference  load  motor,  moreover,  electrical  drive’s  respone 
ensuring  the  “fast,  accurate  and  small  overshoot”  and  reducing  the  resonance 
oscillations  for  two-mass  system  using  induction  motor  fed  by  voltage  source 
inveter  with  ideally  control  performance  of  stator  current.  Backstepping 
controller  uses  the  non-linear  equations  of  an  induction  motor  and  the  linear 
dynamical  equations  of  two-mass  system,  the  Lyapunov  analysis  and  the 
errors  between  the  real  and  the  desired  values.  The  controller  has  been 
implemented  in  both  simulation  and  hardware-in-the-loop  (HIL)  real-time 
experiments  using  Typhoon  HIL  402  system,  when  the  drive  system  operates 
at  a  stable  speed  (rotor  flux  is  constant)  and  greater  than  rated  speed  (field 
weakening  area).  The  simulation  and  HIL  results  presented  the  correctness 
and  effectiveness  of  the  controller  is  proposed;  furthermore,  compared  to  PI 
method  to  see  the  response  of  the  system  clearly. 
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1.  INTRODUCTION 

The  two-mass  system  with  flexible  shaft  can  be  found  in  many  industrial  applications  such  as 
machine  tools,  rolling-mills,  robot  arms,  conveyer  belts,  etc.  Due  to  the  finite  stiffness  and  damping 
coefficients  of  the  shaft,  the  speed  oscillation  of  the  drive  system  is  inevitable  and  may  cause  detrimental 
influences  on  not  only  the  quality  of  the  products  but  also  the  mechanical  and  electrical  components  of  the 
drive  system.  Hence,  speed  control  design  for  these  drive  systems  plays  a  key  role  in  improving  the  system 
characteristic  and  products  quality.  Suppressing  the  shaft  torsional  vibration,  rejecting  the  load  disturbance 
torque  and  tracking  the  reference  speed  are  main  goals  of  the  controller  [1-4].  However,  these  tasks  are  not 
easy  to  be  achieved  due  to  the  fact  that  not  all  system  state  variables  are  available,  furthermore,  unknown 
nonlinearities  such  as  friction  and  backlash  may  also  cause  undesirable  inaccuracy  for  the  control  system. 

Various  control  schemes  have  been  developed  in  the  literature  to  suppress  the  torsional  vibration  of 
the  aforementioned  two-mass  system.  Obviously,  the  classical  cascade  structure  using  proportional-integral 
(PI)  speed  controller  is  the  most  popular  control  scheme.  However,  this  method  cannot  effectively  suppress 
the  vibration  since  the  PI  controller  does  not  have  enough  degree-of-freedom  to  handle  such  high  order 
dynamic  system  like  this  [5].  Thus,  advanced  control  techniques  such  as  sliding  mode  control,  flatness-based 
control,  adaptive  control,  and  fuzzy  adaptive  control  are  also  employed  to  get  a  better  performance  [6-11]. 
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In  this  research,  the  field-oriented  control  (FOC)  which  has  been  widely  used  in  industrial  inverters 
is  used  to  control  the  induction  motor  of  the  drive  system  [12].  Based  on  an  observation  that  if  the  response 
of  the  stator  current  controller  is  fast  enough,  the  order  of  the  two-mass  dynamic  system  can  be  reduced  from 
7th  to  3rd  in  nominal  speed  range  and  from  7th  to  4th  in  field  weakening  range.  Consequently,  the  control 
design  for  the  reduced  order  system  can  be  simplified.  Since  the  dead-beat  controller,  which  forces  the  stator 
current  to  its  desired  value  in  finite  step,  has  already  successfully  been  developed  for  the  stator  current  loop 
in  our  previous  study  [13],  the  remaining  objective  of  this  research  is  to  achieve  a  high  dynamic  speed 
control  with  torsional  oscillation  reduction  for  the  drive  system. 

In  this  paper,  a  backstepping  control  algorithm  is  proposed  to  control  the  load  speed  of  two-mass 
systems.  For  achieving  this  goal,  a  two-block  control  strategy  is  proposed.  It  is  backstepping  control 
algorithm  for  load  speed  and  motor  speed.  On  the  other  hand,  rotor  flux  control  is  designed  acording  to 
backstepping  control  algorithm.  Generating  the  desired  motor  torque  by  dead-beat  current  control,  not  only 
guarantees  accurate  load  speed  but  also  proves  the  stability  of  the  system  [14,  15]. 


2.  MODELING  OF  DRIVE  SYSTEM 

2.1.  Mathematical  model  of  three  phase  induction  motor 

In  field  synchronous  coordinate,  the  three-phase  induction  motor  canbe  described  by  the  following 
dynamic  equations  [12]. 
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In  which,  lgd,lsq 


are  ^  components  of  the  stator  current; 


UJ,UJS 


are  mechanical  rotor  velocity  and 


synchronous  speed,  respectively;  ^ rd  ’  * rq  are  dq  components  of  the  rotor  flux  ;  &  is  total  leakage  factor;  r  is 

rotor  time  constat: Usd  ,  U$q  are  ^  components  of  the  stator  voltage;  is  stator  inductance. 

It  can  be  seen  that  the  original  state  (1)  is  bilinear  and  is  of  4th  order.  Assume  that  the  current 
controller  is  perfect  with  ideal  response,  the  induction  motor  model  can  be  reduced  as: 
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Where: 
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In  which,  m  5  Lr  is  mutual,  rotor  inductance;  lm  is  vector  of  magnetizing  current;  Ml  is  torque  load; 


Z  T 

p  is  number  of  pole  pairs;  J  is  torque  of  inertia; 

The  state  (2)  is  of  2ndorder,  stator  current  isd  is  used  to  control  the  motor  flux  and  isq  is  dedicated  to 
speed  control. 


2.2.  Model  of  the  two-mass  system 

Typical  configuration  of  a  two-mass  system  show  as  Figure  1 : 
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Figure  1 .  Typical  configuration  of  a  two-mass  system 


The  typical  configuration  of  a  two-mass  system  is  illustrated  in  Figure  1.  The  system  can  be 
described  by  the  following  linear  dynamical  equations  [2]. 
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Where: 

’  ^2  Are  motor  speed,  load  speed;  ^  are  motor  and  load  angle  acceleration;  ^  is  angle;  ^  is 
shaft  damping  coefficient;  cis  shaft  stiffness. 

When  the  drive  system  operates  at  field  weakening  range,  the  mathematical  model  of  the 
system  is  [2]. 
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For  control  design  purpose,  term  (4)  is  rewritten  in  the  following  state-space  form: 
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3.  BACKSTEPPING  CONTROL  DESIGN  FOR  THE  FLUX  AND  SPEED  CONTROL 

Based  on  [13]  and  the  concept  of  the  backstepping  design  for  induction  motor  [10],  the  proposed 
control  scheme  for  the  two-mass  system  using  induction  motor  fed  by  VSI  with  perfect  stator  current 
controller  is  shown  in  Figure  2. 


Figure  2.  Block  diagram  control  of  two-mass  system 


In  this  section,  a  backstepping  control  stagegy  is  proposed  to  control  the  rotor  flux  and  speed  load  of 
two-mass  system  with  flexible  couplings.  This  paper  presented  the  induction  motors  operate  at  rated  require 
speed  is  greater  than  the  rated  speed;  therefore,  the  flux  must  be  reduced  to  ensure  that  the  power  supply  is 
not  overloaded. 

3.1.  Flux  controller  design 

Since  rotor  flux  magnitude  are  our  control  variables,  the  tracking  errors  are  defined  by: 

.  .*  .* 

Z]  ~  l,n  l,n  ( lfn  is  the  desired  value  of  flux  )  (7) 

taking  the  derivative  of  both  sides  of  (7)  Gives: 

.  _  di^  J_  J_  di*m 

1  ’  ^dt  Tr  m  Tr  *  dt 


we  can  see  that  (8)  is  the  non-linear  equation 
where: 


dx 

dt 


:/(x,u,0 


z‘  is  the  state  variable;  lsd  is  the  input  signal  (or  control  signal). 

to  stabilize  the  flux  signal  given  by  (8),  the  candidate  Fyapunov  fuction  can  be  defined  as  following: 


r(z,)  =  -tz,2>0 


(9) 


which  the  derivative  of  (9)  is: 
V(z,)  =  z,z,  = 


(10) 
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The  Lyapunov  stability  condition  can  be  verified  if  the  stabilization  function 

follows: 


is  chosen  as 


K(z,)  =  -c1z1J<0Vzl 


(11) 


where  1  is  the  positive  constant  design  that  is  determine  the  closed-loop  dynamics 
Substitute  (11)  into  (10)  Result  in: 


hd  ~  ~C\TrZ\  +  L  +  Tr  ~ 

dt 


with: 


l  zz  i 
sq  *$q 


The  (12)  as  shown: 


=  -ClTrZl+im+Tr^f- 
dt 


(12) 


(13) 


(14) 


so,  the  lsd  is  the  actual  control  signal,  the  control  rule  for  the  flux  controller. 

3.2.  Speed  controller  design 

For  designing  this  nonlinear  control  law,  it  is  assumed  that  all  feedback  signals  from  the  motor,  the 
shaft  and  the  load  sides  are  available  in  this  section.  As  Figure  3  shows,  the  proposed  backstepping  control 
algorithm  is  designed  in  two  steps.  In  each  step  a  subsystem  will  be  studied.  Each  step  includes  the  dynamics 
of  the  previous  subsystems  for  backstepping  technique.  Selecting  a  proper  Lyapunov  function  for  each  step, 
guarantees  the  asymptotic  stability  of  the  related  subsystem.  Completing  all  two  steps,  results  in  not  only  the 
asymptotic  stability  of  the  whole  system,  but  also  the  desired  input  (motor  torque)  law  for  accurate  load 
speed  control. 


Second  step 


Figure  3.  Block  diagram  of  two-step  backstepping  control  algorithm  of  a  two-mass  system 


The  speed  control  equation  is  given: 
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The  (15)  has  triangular  form  as  follows: 


dr 

mo  


df  . .  "  (16) 

these  two  step  for  designing  the  backstepping  control  law  are  mentioned  in  the  following: 

Step  1:  apply  the  backstepping  for  load  speed  *Pi : 

Since  load  speed  magnitude  are  our  control  variables,  the  tracking  errors  are  defined  by: 


Zl  v  u  '  :  C  1  Is  the  reference  load  speed) 
taking  the  derivative  of  both  sides  of  (17)  Gives: 


(17) 


22=^-^  (18) 
to  stabilize  the  flux  signal  given  by  (18),  the  candidate  Lyapunov  fuction  can  be  defined  as  following: 


y(z2)  =  -z2>  0 


which  the  derivative  of  (19)  Is: 

V(22)  =  2Zn 

then: 
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,  1  c  A 

=^j  +  -W4--7A* 
a  a 


(22) 


The  Lyapunov  stability  condition  can  be  verified  if  the  stabilization  function  is  chosen  as 


follows: 


F(24)» 


<0Vr2 


(23) 


where  d  is  the  positive  constant  design  that  is  determine  the  closed-loop  dynamics. 
From  (18)  and  (15)  We  obtain: 

d  i  ...  d  .  c  A  d  .  ,  1  , 
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then 
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so  load  speed  shown  as: 


(24) 


(25) 
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Step  2:  apply  the  hackstepping  for  motor  speed  (^  ]): 

Since  load  speed  magnitude  are  our  control  variables,  the  tracking  errors  are  defined  by: 


"!  >  lcf  1 1  Is  the  reference  motor  speed) 

taking  the  derivative  of  both  sides  of  (27)  Gives: 


(27) 


(28) 


To  stabilize  the  flux  signal  given  by  (28),  the  candidate  Lyapunov  fuction  can  be  defined  as 
following: 


T  /  J_  -  2  _ |_  _j_  2 

1_  2  2  2  1 


which  the  derivative  of  (29)  is: 

Vx  »  Z2  Z2  +  Zj  Zj 

then  substitute  (15)  and  (28)  into  (30),  we  get 
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1  dJx 

i1'  =  T-(cA>p+d<plll-d<P2+—z2) 


(29) 

(30) 

(31) 
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so,  the  sq  is  the  actual  control  signal,  the  control  rule  for  the  speed  controller. 

3.3.  Design  the  reference  trajectories 

When  designing  the  controller,  we  usually  set  the  desired  value  to  a  constant  (step-by-step 
trajectory).  At  the  first  interval  of  time,  the  large  deviation  will  cause  a  large  input  requirement.  However,  the 
response  of  the  actuator  is  finite,  it  is  unlikely  that  the  controller  will  meet  the  requirements.  Thus  the 
backstepping  control  method  will  "soften"  the  desired  signal  by  changing  the  desired  signal  slowly  so  that  the 
response  of  the  system  is  soft  and  better.  For  this  reason,  the  second  order  is  chosen  for  designing  the 
reference  trajectories.  Because  the  control  rules  are  represented  in  (14)  And  (32)  Use  the  first  order 

derivative,  the  form  of  trajectories  l,n  and  u  will  have  to  be  first  order  differentiable. 


Flux  reference  trajectory  designing  l,n 

C  G)  l  l 

idm(s)  ~  (1  +  7>)2  ~  l  +  27>  +  7;Y 

T  i* 

where:  1  is  time  constant;  lyn  is  output  of  flux  reference  trajectory; 


(33) 


m  is  input  of  flux  reference  trajectory. 


Int  J  Pow  Elec  &  Dri  Syst,  Vol.  10,  No.  2,  June  2019  :  720  -  730 


Int  J  Pow  Elec  &  Dri  Syst 


ISSN:  2088-8694 


□  727 


Therefore,  on  the  time  domain: 


d  im 
dt2 


-IT 


dim 

dt 


trajectories  of  flux  reference: 


£L 

dt 7 


(34) 


(35) 


Speed  reference  trajectory  designing  U 
u\s)  _  1  _  1 

ud(s)~  (1  +  T2sf  ~  1  +  2T2s  +  T22s2 


(36) 


where:  is  time  constant;  u  is  output  of  speed  reference  trajectory;  ^  is  input  of  speed  reference 

trajectory.  Therefore,  on  the  time  domain: 


d2u 

dt2 


+  27; 


duo 

dt 


-  +  cj  =  va 


trajectories  of  flux  reference: 


(37) 


(38) 


4.  RESULT  AND  ANALYSIS 

In  this  section,  the  proposed  control  strategy  is  verified  by  simulation  using  Hardware-in-the-Loop 
(HIL).  The  proposed  control  is  also  verified  in  Hardware-in-the-Loop  (HIL)  environment  using  a  Typhoon 
device  in  Figure  4.  This  device  consists  of  an  HIL402  card  that  simulates  grid  source,  load,  and  three  full- 
bridge  with  a  common  DC-link  capacitor  using  IGBTs.  The  system  hardware  is  simulated  in  real  time  on  the 
HIL  platform  with  a  time  step  of  1  ps  that  means  very  nearly  physical  model,  while  the  pulse  width 
modulation  (PWM)  carrier  frequency  is  5  kHz.  Voltage  and  current  controllers  as  well  as  PLL  are 
implemented  in  DSP  TMS320F2808  card. 

The  simulation  parameters  for  rated  power:  0.5kW;  rated  phase  voltage:  220  VRMS;  rated 
frequency:  50Hz;  d=0.313Nm/rad/s;  c=300000Nm/rad. 

Simulation  procedure  is  as  follows: 

At  t  =  Os,  the  magnetic  current  is  created.  Then  at  time  instance  t  =  0.4s,  the  motor  starts  to  speed  up 
to  2500rpm.  At  t=l s,  the  full  load  is  inserted.  At  t  =  1,4s,  the  motor  is  continued  to  speed  up  to  3500rpm  with 
full  load.  Finally,  the  motor  starts  to  change  the  rotating  direction  with  a  reference  speed,  i.e,  2500rpm,  at 
time  instance  t  =  2.0s. 

The  simulation  results  of  the  proposed  method  are  presented  and  compared  to  those  getting  from  the 
conventional  PI  controller  to  evaluate  the  advantage  of  the  backstepping  method. 

The  responses  of  magnetizing  current  are  shown  in  Figure  5.  Both  conventional  PI  and  backstepping 
methods  show  stable  operation  not  only  in  nominal  but  also  field  weakening  range.  However,  it  can  be 
observed  that  the  backstepping  gives  better  transient  response,  i.e.,  without  overshoot  as  seen  in  Figure  5  and 
shorter  settling  time  (0.2s  in  comparison  with  0.35s  for  PI  control). 

The  response  of  the  torque,  speed  as  well  as  torque  ripple  of  the  motor  are  shown  in  Figure  6.  As 
can  be  observed,  the  torque  and  speed  quickly  match  their  commands  for  both  methods.  Nevertheless,  the 
backstepping  always  shows  better  performance,  i.e,  lower  torque  ripple  and  shorter  accelerating  time,  as  can 
be  seen  in  Figure  7. 
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Speed  control  method 


Flux  response 


Speed  control  method 


Speed  response 


t  orque  response 


Figure  4.  Picture  of  hardware  HIL  platform. 
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Figure  6.  Speed  and  torque  responses 
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Figure  7.  Chart  of  torque  and  speed  responses,  (a)  torque  responses,  (b)  speed  responses 


5.  CONCLUSION 

In  this  paper,  the  two-mass  system  with  flexible  couplings  comprises  of  an  induction  motor  and  a 
load  is  considered.  At  the  motor  side,  the  current  response  is  assumed  to  be  ideal  leading  to  a  simplified 
model.  Backstepping  control  is  employed  to  solve  flux  and  speed  control  problem  of  the  system.  The  system 
can  operate  at  field  weakening  region.  The  simulation  results  show  that  high  dynamic  and  suppression  of  the 
mechanical  oscillation  of  the  drive  system  can  be  achieved. 
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